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ABSTRACT 

We report on the detection of a planetary companion in orbit around the primary star of the binary 
system 7 Cephei. High precision radial velocity measurements using 4 independent data sets spanning 
the time interval 1981-2002 reveal long-lived residual radial velocity variations superimposed on the 
binary orbit that are coherent in phase and amplitude with a period or 2.48 years (906 days) and a semi- 
amplitude of 27.5 ms _1 . We performed a careful analysis of our Ca II H & K S-index measurements, 
spectral line bisectors, and Hipparcos photometry. We found no significant variations in these quantities 
with the 906-d period. We also re-analyzed the Ca II A8662 A measurements of Walker et al. (1992) which 
showed possible periodic variations with the "planet" period when first published. This analysis shows 
that periodic Ca II equivalent width variations were only present during 1986.5 - 1992 and absent during 
1981-1986.5. Furthermore, a refined period for the Ca II A8662A variations is 2.14 yrs, significantly 
less than residual radial velocity period. The most likely explanation of the residual radial velocity 
variations is a planetary mass companion with M sin i — 1.7 Mj up it er and an orbital semi-major axis of 
02 = 2.13 AU. This supports the planet hypothesis for the residual radial velocity variations for 7 Cep 
first suggested by Walker et al. (1992). With an estimated binary orbital period of 57 years 7 Cep is the 
shortest period binary system in which an extrasolar planet has been found. This system may provide 
insights into the relationship between planetary and binary star formation. 

Subject headings: planetary systems — techniques: radial velocities 



1. INTRODUCTION 

The first high-precision radial velocity survey for plan- 
etary companions to nearby stars was conducted with the 
use of a HF gas absorption cell on the Canada-France- 
Hawaii Telescope (Campbell & Walker 1979; Walker et al. 
1995). The 7 Cephei system (HR 8974 = HD 222404 = 
HIP 116727) was one of 16 stars on the observing list of 
this program. Campbell, Walker, & Yang (1988) reported 
that two of the stars on the observing list (7 Cephei and 
X 1 Orionis) were previously unknown single-lined spectro- 
scopic binaries. In the case of 7 Cephei, Campbell, Walker, 
& Yang (1988) found evidence for radial velocity "bumps" 
superimposed on the large amplitude binary motion. The 
residuals after removing the binary orbit yielded a period 
of about 2.7 years and a velocity semiamplitude of about 
25 ms -1 . The authors examined several possible causes 
of the observed radial velocity variability, and finally con- 
cluded that the system had a "probable third body." If 



substantiated, this would have been the first detection of 
an extrasolar planetary system. However, Bohlender, Ir- 
win, & Yang (1992) classified 7 Cep as a KO III star, rais- 
ing the possibility that the observed radial velocity vari- 
ations were due to the recently discovered (Walker et al. 
1989; Hatzes & Cochran 1993) long-period radial veloc- 
ity variability of most K giants. A detailed analysis of 
the CFHT radial velocity data on 7 Cep by Walker et 
al. (1992) showed a clear low-amplitude (27ms _1 ) sig- 
nal with a 2.52 year period superimposed on the binary 
orbital motion of indeterminate period. They also found 
evidence of a possible variation of the Ca II A8662 A emis- 
sion line index with the same 2.5 year period, leading them 
to conclude that the observed low-amplitude radial veloc- 
ity variation was most likely due to K-giant variability at 
the period of the stellar rotation. 

Here we present new high-precision radial velocity data 
on the 7 Cephei system obtained from McDonald Obser- 
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vatory. When combined with the CFHT data, we show 
that the 2.5-year low- amplitude radial velocity variability 
of 7 Cep has remained constant for over 20 years. A si- 
multaneous orbital solution for both the binary star and 
the 2.5-year low amplitude variability is computed. We 
demonstrate that there is no correlation between the low- 
amplitude radial velocity variations and the Mt. Wilson 
Ca II S-index, nor is there any indication of photospheric 
absorption line profile variability. These results, combined 
with the current best classification of 7 Cep as a Kl IV star 
(Fuhrmann 2003), all indicate that the preferred interpre- 
tation of the data is that 7 Cep A has a planetary-mass 
companion in a 2.5 year period orbit. 

2. THE RADIAL VELOCITY DATA SETS 

We consider four independent sets of high precision ra- 
dial velocity data for 7 Cephei, covering the interval from 
1981 through 2002. The first of these is from the CFHT 
survey Campbell, Walker, & Yang (1988), Walker et al. 
(1995), with the velocities taken from Table 1 of Walker 
et al. (1992). All of the other data were from the Mc- 
Donald Observatory Planetary Search (MOPS) Program 
(Cochran & Hatzes 2000). Phase I of the MOPS program 
used the telluric O2 lines near 6300A as the velocity met- 
ric, a technique suggested by Griffin & Griffin (1973). A 
single order of the 2.7m coude spectrometer 6-foot camera 
with the echelle grating was isolated onto a Texas Instru- 
ments (TI) 800 x 800 CCD at R = 210, 000. This system 
gave 15-20 ms _1 precision on stars down to about V = 6, 
but suffered from systematic velocity errors, most likely 
due to prevailing atmospheric winds. In 1992 the pro- 
gram switched to a temperature stabilized I2 cell (Koch 
& Wohl 1984; Libbrecht 1988) as the velocity metric for 
Phase II of the MOPS (R = 210, 000 for this data as well). 
This eliminated the systematic errors, and gave a routine 
radial velocity precision of 15 ms -1 . This precision was 
limited by the 9.6 A bandpass of the spectrum, and by 
the poor charge-transfer and readout properties of the TI 
CCD. To solve these problems, and to achieve substan- 
tially improved precision, we began Phase III of the ra- 
dial velocity program in July 1998, using the I2 cell with 
the newly installed 2dcoude cross-dispersed echelle spec- 
trograph (Tull et al. 1994). This instrument when used 
with a Tektronix 2048x2048 detector provides a nominal 
wavelength coverage of 3600 A - 1/zm at a resolving power 
of R = 60, 000. The Tektronix CCD also had significantly 
better charge transfer and readout properties than the TI 
device. The complete spectral coverage of 2dcoude gives 
us two advantages: first, we can utilize the full reference 
spectrum of the l2-cell for the RV determination and sec- 
ond, we can simultaneously determine the stellar chromo- 
spheric emission in the cores of the Call H&K lines to use 
as stellar chromospheric activity indicators. We will dis- 
cuss in detail the question whether there is a correlation 
between the RV-results and the activity indices for 7 Cep 
in section §4. 

To extract the RV-information from the I2 self- 
calibrated spectra taken during Phase III we employed 
our Austral RV-code, which uses a Maximum Entropy 
Method deconvolution in order to obtain a higher resolved 
stellar template spectrum and several reconstruction al- 
gorithms which model the shape and symmetry of the 



instrumental profile (IP) at the time of observation. A 
detailed description of the Austral code can be found in 
Endl, Kiirstcr, & Els (2000). The algorithm follows in 
general the modeling idea first outlined by Butler et. al. 
(1996), and the IP reconstruction techniques by Valenti et 
al. (1995). Using the 2dcoude spectrometer in I 2 -cell self- 
calibration mode and the Austral code for the analysis, 
we obtain a long term RV precision of 5 — 15 ms _1 on a 
routine basis for stars down to a magnitude of V = 9.0. 

All of the MOPS data from Phases I, II and III are given 
in Tables 1-3. The uncertainties quoted there for Phase 
II and III data are the "internal" errors, as represented by 
the rms of the individual spectral chunks about the mean 
value. We regard these as a lower limit on the actual un- 
certainties, since these values do not include the effects of 
any residual systematic errors that may be present. It is 
difficult to estimate internal errors for the Phase I data be- 
cause these are not analyzed in "chunks" like the Phase II 
and III measurements. Phase I observations of a constant 
star (r Cet) show an rms scatter of 23 ms _1 . The rms 
scatter of the Phase I measurements about the final orbital 
solution is 17.4 ms -1 . We list a slightly higher value of 19 
ms as the "error" of the Phase I measurements in Ta- 
ble 1. This error was the most reasonable error to assign 
based upon \ 2 tests of the data and was the value assumed 
for the Phase I measurements in the final orbital solutions. 
The true error for an individual Phase I measurement is 
almost surely higher due to different signal-to-noise ratios 
of individual spectra and systematic errors due to wind 
and temperature and pressure changes in the earth's at- 
mosphere. 

The RV measurements for all data sets are shown in 
Figure 1. A different velocity offset (see below) had to be 
applied to each data set so that they would all have the 
same zero point. One can clearly see that the "wiggles" 
superimposed on the binary variations, that were first re- 
ported by Walker et al. (1992) are still present in the final 
data set (Phase III) taken 20 years later. 

Figure 2 shows a sequence of Lomb-Scargle peri- 
odograms (Lomb 1976, Scargle 1982) of the RV measure- 
ments after subtraction of the velocity variations due to 
the binary companion (see below). The top panel is only 
for the CFHT data. The central panel is for the CFHT + 
McDonald data, excluding the Phase III data. The lower 
panel is the periodogram for all of the RV measurements. 
The increase in power at a period of « 2.5 yrs with the 
addition of each data set is the first indication that these 
residual RV variations are long-lived and coherent. (The 
false alarm probability of the peak for the full dataset is 
» 10- 20 ). 

3. ORBITAL SOLUTIONS 

The program GaussFit (Jcffcrys et al. 1988; McArthur 
et al. 1994) was used to fit simultaneously all orbital pa- 
rameters, both for the stellar and the presumed sub-stellar 
companion, using nonlinear least squares with robust esti- 
mation. Because each instrument produces relative radial 
velocities with their own arbitrary zero point, the indi- 
vidual velocity offsets were a free parameter in the least 
squares solution. Table 4 gives the velocity offsets of the 
different data sets. Subtracting these velocity offsets from 
an individual data sets will place them all on the same 
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radial velocity scale. The line in Figure 1 shows the com- 
bined orbital solution to 7 Cep. Figure 3 shows the "planet 
only" orbital solution after subtracting the contribution of 
the binary orbit to the RV measurements. 

Table 5 lists the orbital parameters for the the stellar 
companion. The errors listed are the correlated errors 
that are produced by GaussFit. The largest error occurs 
in the period since our observations span about 20 years, 
or about one-third of the binary orbit. Note that the re- 
duced x 2 °f t ne binary "only" orbital solution is rather 
large (4.36) indicating the presence of additional RV vari- 
ations. 

Griffin et al. (2002) presented an orbital solution for the 
stellar companion to 7 Cep. They combined radial velocity 
measurements spanning over 100 years taken at five obser- 
vatories. These measurements included the CFHT data as 
well as some of our McDonald Observatory measurements 
(read from an enlarged copy of a published graph!). The 
parameters of their orbital solution is also listed in Table 5. 
In spite of the more limited time span of our observations 
our orbital solution agrees quite well with the Griffin et al. 
orbit. 

Table 6 lists the GaussFit orbital parameters for the 
sub-stellar companion. Also listed is the rms scatter of 
the individual datasets about the combined orbital solu- 
tion. The period of 2.48 years and semi-amplitude of 27.5 
ms _1 are consistent with the values found by Walker et 
al. (1992) using only the CFHT data set. Note that the 
reduced x 2 is significantly lower ( x 2 = 1-47) when includ- 
ing the planet in the orbital solution. Fuhrmann (2003) 
estimates the primary mass of 7 Cep as M = 1.59 ± 0.12 
Mq. This results in a planetary mass of M p sini = 1.7 ± 
0.4 M j U piter and an orbital semimajor axis of 2.13 AU. 

Figure 4 shows the phase diagram of the individual data 
sets. The phase, amplitude, and overall shape of the RV 
curves hold up remarkably well over a 20 year time span, 
again evidence for long-lived and coherent variations that 
are consistent with the presence of a planetary companion. 

4. THE NATURE OF THE RV VARIATIONS 

The question naturally arises as to whether some phe- 
nomena instrinsic to the star (spots, convective shifts, pul- 
sations) are responsible for the 2.48 yr period. After all, 
7 Cep is most likely a subgiant (see below), a class of stars 
for which magnetic activity, surface structure, convection, 
and pulsations are poorly known. Moreover, Walker et 
al. (1992) did find evidence for weak periodic variations 
in the equivalent width of the Ca II A8662 A line with the 
same period as the planet. In spite of the intrinsic variabil- 
ity of giant stars planetary companions have been found 
around other giant stars (Frink et al 2002; Setiawan et al. 
2003) and one subgiant (Butler et al. 2001). So planets 
around giant evolved stars is not an unreasonable hypoth- 
esis. Before concluding that the planet hypothesis is the 
most likely explanation for the residual RV variations we 
must demonstrate that 7 Cep does not exhibit any signifi- 
cant variations with the planet period in other quantities. 
Here we examine whether 7 Cep also has spectral and pho- 
tometric variations. 

4.1. Bisector Analysis 



The spectral line shapes can also provide evidence in 
support of the planetary hypothesis. Surface features 
(Hatzcs 2002) or nonradial pulsations (Hatzes 1996, Brown 
et al. 1998) should produce changes in the spectral line 
shapes with the same period as the RV variations. A con- 
venient means of measuring the asymmetry of a spectral 
line is the line bisector, or the line segments connecting 
the midpoints of the spectral line from the core to the 
continuum. Line bisector studies have been used to es- 
tablish the planetary nature of 51 Peg (Hatzes, Cochran, 
& Johns-Krull 1997; Hatzes, Cochran, & Bakker 1998a,b) 
and the starspot nature of the planet-like RV variations in 
HD 166435 (Queloz et al. 2001). 

The McDonald Observatory Phase I observations taken 
in the 6300 A region provide an excellent data set for 
studying possible line profile variations in 7 Cep. The 
data has very high spectral resolution (R = 210,000) and 
there are several strong spectral lines that were velocity 
shifted clear of the telluric features because of the Earth's 
barycentric velocity. The Fe I A6301.5 A feature was free 
of telluric lines over the entire data set, Fe I A6297.8 for 30 
observations, and Fe I A6302.5 A for 8 observations. Both 
the bisector velocity span and curvature were measured 
for these lines. The velocity span is defined as the velocity 
difference between two arbitrary points on the line bisector 
while the curvature is the difference between the velocity 
span of the top half of the bisector minus the velocity span 
of the lower half of the bisector. 

A 5th order polynomial was fit to each measured line bi- 
sector and the velocity span points taken at 0.4 and 0.8 of 
the continuum. (In measuring the line bisector one should 
avoid both the core and the continuum where the errors in 
the bisector measurement become large.) The additional 
point required for the curvature measurements was taken 
at 0.6 of the continuum value. The mean value of the bisec- 
tor span and curvature for each line was then subtracted 
from the individual measurements and the residual span 
and velocity measurements were then averaged, weighted 
by the rms scatter of the bisector measurements for each 
spectral line. 

Figure 5 shows the phase-binned averages (A(f> « 0.05) 
of the bisector span and curvature variations for 7 Cep 
phased to the planet orbital period. The short-dashed 
horizontal line gives the zero-point reference and the long- 
dashed lines show the velocity extrema of the residual RV 
variations due to the planetary companion. There are no 
obvious phase variations of the bisector quantities with 
the planet period. The least squares sine fit to the bisector 
quantities assuming a period of 906 days yields amplitudes 
of 4.8 ± 4.4 ms _1 and 2.4 ± 4.5 ms _1 for the bisector and 
span variations, respectively. These amplitudes are consis- 
tent and are significantly less than the observed 27 ms _1 
radial velocity variations of the planet. 

The McDonald Phase I data had the largest errors of all 
the data sets and the phase variations of these show the 
least evidence for the 906 day period (Fig. 4). One could 
argue that there are no residual RV variations with which 
to correlate to the bisector measurements. However, the 
Phase I data is still consistent with the presence of the 
planet signal. Not only is this signal present in the Phase 
I and II data alone, but the power increases by almost a 
factor of two over the individual periodograms when the 
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two data sets are combined. Although not obvious to the 
eye, the 906-day period is still present in the Phase I data. 
On the other hand, the averaged, phase-binned bisector 
measurements have variations significantly less than the 
906-day RV amplitude and we believe this adds additional 
evidence (along with the photometric and Ca II analysis 
presented below) in support of the planet hypothesis. If 
the bisector variations had scatter comparable to the RV 
amplitude, then one could at least make a plausible ar- 
gument in favor of possible bisector variability. Figure 5 
excludes that. 

4.2. Photometric Variations 

The Hipparcos mission (Perryman 1997) took high pre- 
cision photometric measurements of 7 Cep contemporane- 
ously with the measurements used for our RV study. If 
cool spots on the stellar surface, or some form of stellar 
pulsations were causing the residual RV variations, then 
this should be evident in the Hipparcos photometry 

Figure 6 shows the Hipparcos photometry taken from 
1989.9 to 1993.2 phased to the planet orbital period. 
Crosses represent individual measurements while solid 
points represent phase-binned averages (A<j> < 0.05). Error 
bars represent the rms scatter of the measurements used 
for each bin. There arc no obvious photometric variations 
with the planet period as confirmed by a periodogram 
analysis of the daily averages of the Hipparcos photometry 
(Figure 7). A least squares sine fit to the phased photo- 
metric data yields an amplitude of AV" = 0.001 ± 0.0009 
mag for any possible photometric variations. (A fit to the 
phase-binned measurements yields AV = 0.0003.) The 
lack of photometric variations in 7 Cep is also consistent 
with the planet hypothesis for the residual RV variations 
for this star. 

4.3. Ca II Variations 
4.3.1. S-index Measurements 

Stellar activity in 7 Cep could induce significant peri- 
odic centroid shifts in photospheric absorption lines which 
could be confused for perturbations made by planetary 
companions (Saar & Donahue 1997; Saar & Fischer 2000, 
Qucloz et al. 2001). For the Phase III data an instrumen- 
tal setup was chosen so as to include the Ca 11 H and K 
lines on the detector. 

To measure stellar chromospheric activity, the Mt. Wil- 
son S index was adopted. This index is defined (e.g. Bal- 
iunas et al. 1992) as a quantity proportional to the sum 
of the flux in 1 A FWHM triangular bandpasses centered 
on the Ca II H and K lines divided by the sum of the flux 
in 20A bandpasses in the continuum at 3901 and 400lA 
(Sodcrblom, Duncan, & Johnson 1991). 

The four quantities to be measured (the two calcium 
line core fluxes plus the two continuum bandpass fluxes) 
are spread across 3 echelle spectral orders which overlap 
by several A. In order to be consistent with activity mon- 
itoring with our other programs (Paulson et al. 2002), 
we did not use measurements in the bluemost order, i.e. 
the continuum region centered on 390lA. In addition, we 
did not measure the Ca II H line (at 3968. 47A) because 
the wings of strong Balmer He feature (at 3970. 07A) are 
within the measured Ca 11 bandpass. This can adversely 



affect the measurement of the Ca 11 H line flux. Therefore, 
we have defined an index SmcD which is the ratio of the 
flux in a 1 A triangular bandpass centered on the Ca II K 
line to the flux in a 20 A bandpass centered on the red- 
ward continuum at 4001 A. Thirty of our program stars 
(for the McDonald Observatory Planet Search) have pre- 
viously been measured as part of the Mt. Wilson survey 
(Baliunas et al. 1995; Duncan et al. 1991). We find a 
linear relationship between our measurements and those 
of the Mt. Wilson survey of the form: 

^Mt.wuson = 0.038(±0.006) + 1.069(±0.040) x S Mc u (1) 

Both the slope and the intercept of the formal fit are within 

I o~ of SmcD = ^Mt. Wilson- We are thus able to transform 
our data into a standard Mt. Wilson S index scale using 
Equation (1). 

The periodogram of the S-index measurements (Fig- 
ure 8) shows no significant power at the orbital frequency 
of the planet. The false alarm probability of the highest 
peak, assessed using a bootstrap randomization process 
(Murdoch et al. 1993; Kiirster et al. 1997), is 50%. The 
lack of variability at the planet period is substantiated by 
phasing these variations to the 906 day period (Figure 9). 
Clearly no significant sinusoidal variations in this chromo- 
spheric index are present. The transformed Ca II K S- 
indcx for 7 Cep is —5.3 which implies a level of magnetic 
activity less than that of the Sun. 

4.3.2. Analysis of the Ca II \8662 A data from Walker et 

al. (1992) 

A conclusion regarding the planetary nature of the resid- 
ual RV variations of 7 Cep would not be complete with- 
out a discussion of the Ca II variations found by Walker 
et al (1992). Although the McDonald S-index measure- 
ments do not show evidence for rotational modulation, Ca 

II A8662 A equivalent width measurements (W\) of Walker 
et al. (1992) did show a hint of sinusoidal variations when 
phased to the planetary orbital period. Although this 
on its own does not completely refute the existence of a 
planet, having Ca II variations with the same period would 
cast more doubt on this hypothesis. For these reasons we 
made a careful examination the significance of the Ca II 
variations reported by Walker et al. (1992). 

Figure 10 shows the correlation between the RV and the 
changes in the Ca II A8662A equivalent width. The two 
quantities show no obvious correlations. The correlation 
coefficient, r, is only 0.08 and the probability that the two 
quantities are not correlated is 0.52. However, there are 3 
obvious outliers in the figure (|AWa| > 3%). Eliminating 
these 3 points increases the correlation coefficient to r = 
0.23 with a probability of 0.08 that the quantities arc un- 
correlated. However, this is still not a strong correlation. 

The McDonald S-index data also do not seem to be cor- 
related to the RV measurements. Figure 11 shows the 
correlation between the S-index and radial velocity mea- 
surements from the McDonald data. The correlation co- 
efficient in this case is r = —0.3 and the probability that 
the two quantities are uncorrelated is 0.13. 

Figure 12 shows the Lomb-Scargle periodogram of the 
variations in the Ca II A8662 A equivalent width measure- 
ments of Walker et al. (1992). (There was a clear outlier in 
the Walker et al. Ca II data with AEW = -4.3. This was 
also evident in Fig. 3 of Walker et al. (1992). This data 



Planet around 7 Cephei A 



5 



point was eliminated before performing the periodogram 
analysis.) Although the most power occurs at a periods 
of w 15 days (top panel), there appears to be significant 
power at the planet period in the expanded scale in the 
lower panel (the vertical line indicates the location of the 
planet RV frequency). This was the basis for the Walker 
et al. favoring rotational modulation as the cause of the 
RV variations. 

A more detailed examination of the CFHT Ca II 
A8662 A equivalent measurements show that the variation 
in these are not long-lived and may have no relation to the 
observed 906-day RV period. We divided the CFHT Ca 
II data into two data sets, one spanning 1981-1986.5 and 
the second 1986.5-1992. The division of the data set was 
taken so as to maximize the power in the long period vari- 
ations found in the second set. (Our results do not change 
substantially if we divide the data set into two equal num- 
bered points.) This resulted in 21 data points in the first 
data set and 29 in the second set. 

Figure 13 shows the periodograms of these two data sets. 
The 1981-1986.5 set shows no power in the frequency in- 
terval < v < 0.01 cd" 1 . The 1986.5-1992 data set shows 
significant power, but at a frequency corresponding to a 
period of 781 ± 116 days, significantly less than the planet 
orbital period. Extending the periodogram to higher fre- 
quencies shows that this is the highest peak out to the 
Nyquist frequency, in contrast to the case for the full data 
set where the low frequency feature was the second high- 
est peak. A periodogram analysis of the full Ca II data 
set lowers the Lomb-Scargle power near v = 0.011 cd _1 
from about 9.0 to 6.7 and increases the period slightly. 
Unlike the case for the RV data, increasing the number 
of measurements does not increase the power in the peri- 
odogram of the feature of interest. Clearly this signal is 
not long-lived. 

The statistical significance of this signal was examined 
using the bootstrap randomization technique. The Ca 
II A8662 A equivalent width measurements over the time 
span 1986.5 - 1992 were randomly shuffled keeping the ob- 
served times fixed. A periodogram was then computed for 
each "random" data set. The fraction of the periodograms 
having power higher than the data periodogram in the 
range 0.0005 < v < 0.01 cd _1 is the false alarm probabil- 
ity that noise would create the detected signal. After 10 5 
shuffles there was no instance of a random periodogram 
having power higher than that found in the data peri- 
odogram. The false alarm probability is thus < 10~ 5 . This 
signal, in spite of only being present in the data for the 
last half of the data set, is highly significant. 

The period of the Ca II A8662A measurements is sig- 
nificantly less, by la, from the residual RV period. Fig- 
ure 14 shows the phase diagrams of the Ca II A8662A 
measurements. The top panel are the measurements dur- 
ing 1986.5-1992 phased to the 781 day period found in the 
Ca II data. The central panel is the same data phased to 
the 906 day residual RV period. Although there are slight 
sinusoidal variations when phased to the 906 day period, 
phasing to the 781 day period produces significantly less 
scatter. By comparison, the lower panel shows the Ca II 
measurements from 1981 - 1986.5 phased to the 906 day 
period. Crosses are individual measurements whereas solid 
points represent phase-binned averages (the error bar rep- 



resents the scatter of the data used for each phase bin). 
There are no significant sinusoidal variations in the 1981 

- 1986.5 data set. 

A least squares sine fit to the the phased Ca II yields an 
amplitude of 1.67 ± 0.25 mA for the dates spanning 1986.5 

- 1992 (1.03 ± 0.25 assuming a 906 day period). For the 
dates covering 1981 - 1986.5 the amplitude of the £Ws662 
variations is 0.48 ± 0.33 mA. Clearly there are amplitude 
variations in the Ca II A8662. 

We conclude that the periodic variations in Ca II found 
by Walker et al. (1992) arc not long-lived and were only 
present during 1986.5-1992. This period was clearly not 
present in our Ca II S-index measurements spanning 1998- 
2002. Since the RV variations are unchanged during the 
entire time span 1981-2002 we do not believe that the 906 
day RV period is related to the cause of the Ca II A8662A 
variability. 

4.4. The Spectral Classification of 7 Cep 

One reason that Walker et al. (1992) favored rotational 
modulation for the residual RV variations was the reclas- 
sification of 7 Cep as a K0III. This was based on a visual 
comparison of the spectrum of 7 Cep to other K giants 
(Bohlender et al. 1992). Recently, Fuhrmann (2003) pre- 
sented an extensive spectral analysis of nearby stars in the 
galactic disk and halo. Included in this study was 7 Cep. 
Table 7 list the stellar properties derived by Fuhrmann. 
Also listed is the Hipparcos distance {dm P ) and distance 
determined spectroscopically (ds P ) using the derived stel- 
lar parameters. These distances agree to within 3% indi- 
cating that we can have some confidence in the results of 
the spectral analysis. Fuhrmann's classification of K1IV 
for 7 Cep is consistent with the sub-giant status for 77 
Cep (K0IV) which has a comparable effective tempera- 
ture (4990 K), gravity (log g = 3.4), radius (4.14 i? Q ), and 
bolometric magnitude (M^oi — 2.3). The most current and 
best analysis of 7 Cep supports this being a subgiant star. 

5. SEARCHING FOR ADDITIONAL COMPANIONS 

The long time-baseline of RV-monitoring of 7 Cep and 
the emergence of several multiple extrasolar planetary sys- 
tems from Doppler surveys (e.g. Butler et al. 1999), en- 
couraged us to search for additional companions in this 
system. For this purpose we performed a period search 
within the RV-residuals, after subtracting both the binary 
orbital motion and the first planetary signal. 

Fig. 15 shows the Lomb-Scargle periodogram of the RV- 
residuals in the period range of 2 to 7900 days (frequency 
= 0.00013 - 0.5 cd _1 ). The strongest peak is found at 
P» 11 days but its significance level is very low (false- 
alarm-probability (FAP) is 6.5%). The FAP-levels shown 
in the figure were determined by 10, 000 runs of a boot- 
strap randomization scheme. 

It is clear from this analysis that no additional periodic 
signal above the noise level is present in the RV-residuals 
of 7 Cep and we thus conclude that the 906-day period 
planet is the only giant planet in this system evident in 
our data.. 

6. DISCUSSION 

Precise stellar radial velocity measurements of 7 Cep 
from 4 independent data sets spanning over 20 years show 
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long-lived, low-amplitude RV variations (K = 27.5 ms -1 ) 
superimposed on the larger radial velocity variations due 
to the reflex motion caused by a stellar companion. We in- 
terpret the low amplitude, shorter period variations as due 
to the presence of a planetary companion with M sin i = 
1.7 Mj up it er and an orbital semi-major axis, a = 2.13 AU. 
Our conclusion that these short period variations are not 
due to rotation, pulsations, or changes in the convection 
pattern of the star is based on several facts: 

1. The 2.48 yr period has been present for over 20 
years with no changes in phase or amplitude during 
this time. 

2. No variations with this period are seen in the Mc- 
Donald Ca II S-indcx measurements. 

3. Spectral line bisector span and curvature measure- 
ments for 7 Cep that are constant to less than 5 
ms _1 over an orbital cycle of the planet. 

4. Contemporaneous Hipparcos photometry that is 
constant to less than 0.001 mag over an orbital cy- 
cle. 

5. The periodic variations in the Ca II A8662 A found 
by Walker et al. (1992) were only present during 
1986.5-1992 and are thus most likely not associ- 
ated with the residual RV variability observed for 
this star. 

The CFHT Ca II data do show evidence for long pe- 
riod variability. Although this signal is weakly detected 
in the full data set, it is much stronger in only the last 
half of the data set spanning 1986.5-1992, and completely 
absent in the data from 1981 - 1986.5. Furthermore, our 
S-index measurements made during 1997 - 2003 show no 
modulation. If rotational modulation was responsible for 
the residual RV variations then we would have seen RV 
amplitude changes. However, the amplitude and phase for 
the residual RV variations remained constant. It therefore 
seems unlikely that rotational modulation is the cause of 
the 906 day period. Furthermore, the best fit period to 
the long term Ca II variations during 1986.5 - 1992 is 781 
± 116 days, which is significantly less than the 906 day 
RV period. 

Gamma Cep shows no Ca II modulation during epochs 
centered on 1994 and 2000, and modulation (781 day pe- 
riod) during a time centered on 1989. This indicates a 
possible "activity cycle" period of 10 - 15 years, much 
longer than the presumed planet period. 

Our spectral line bisector measurements, which are con- 
stant, also support our conclusion that we are not seeing 
rotational modulation. One could argue that the residual 
RV variations are due to changes in the convection pat- 
tern in the star and not due to magnetic structure (plage, 
spots). For instance, if the ratio of areas of convective, 
hot rising cells and intergranule, sinking lanes changes 
with an activity cycle (in this case with a period of ~ 
900 days), then the amount of convective blue (red) shift 
would change periodically resulting in a measured RV sig- 
nal. If the stellar magnetic fields are strong enough to alter 
the convection pattern of the star, but too weak to cause 
chromospheric structure, then we would see RV variations 



without strong variations in the Ca II emission. However, 
in this case we should also see changes in the spectral line 
bisectors. The lack of significant variability in the spectral 
line bisectors does not support the hypothesis of a chang- 
ing convection pattern on the star, or a least changes that 
can influence the RV measurements. 

The inclination angle of the rotation axis can be esti- 
mated by comparing the expected equatorial rotational ve- 
locity to the projected rotational velocity (vsini), assum- 
ing that 781 days (period of the Ca II A8662 A variations) 
is the true rotation period of 7 Cep. Fuhrmann (2003) 
estimates a radius R — 4.66 Rq which gives an equato- 
rial rotational velocity of 4.9 kms -1 . We have measured a 
vsini = 1.5 ± 1.0 kms -1 which is consistent to the value 
determined by Fuhrmann (2003). This yields a sini = 0.5 
-0.1. Assuming that the orbital and rotation axis of the 
star is aligned, the true mass for the planetary companion 
to 7 Cep is w 3-16 M Jupiter . 

Given that the 906-day RV variations are due to a sub- 
stellar companion, the question arises whether this system 
is stable. After all, both the planet and the stellar com- 
panion have modestly eccentric orbits. Preliminary indi- 
cations are that this system is indeed stable (Dvorak et 
al. 2003). (The study by Dvorak et al. was made prior 
to the publication of our paper and it used preliminary 
orbital parameters published only in conference abstracts. 
The stability analysis should be re-done using the final 
parameters and error estimates we have presented here.) 

The 7 Cep system presents a very interesting system for 
the study of planet formation. Planets have been found 
around host stars that are in a binary system, but these 
are widely separated pairs so the presence of the stellar 
companion may not have an influence on the planet for- 
mation around one of the stars. The 7 Cep binary is the 
shortest period binary for which a planetary companion 
has been found in orbit around one of the components. 
This implies that even in such a relatively close binary 
the presence of the stellar companion does not hinder the 
process of planet formation. 

One final comment. RV surveys have had a stunning 
success at finding planets in orbit around other stars. Yet 
this is still an indirect detection method with the disad- 
vantage that various stellar phenomenon can mimic a plan- 
etary signal. RV searches for extrasolar planets are, and 
should continue to be a careful process with many candi- 
dates undergoing the painstaking process of confirmation. 
Several extrasolar planet candidates have a proved to be 
due to stellar phenomena (Qucloz et al. 2001; Butler et al. 
2002) and more false detections will undoubtedly be un- 
covered in the future. Likewise, radial velocity variations 
that were initially attributed to stellar phenomena in the 
past may ultimately prove to be due to a planetary com- 
panion, as is apparently the case for 7 Cep. RV searches for 
extrasolar planets must be accompanied by careful studies 
of the host star and the measurement of a various pho- 
tometric and spectroscopic (Ca II emission, line shapes) 
quantities. These measurements are important for they 
could have confirmed the planet around 7 Cep seven years 
prior to the discovery of 51 Peg b. 

This material is based upon work supported by the Na- 
tional Aeronautics and Space Administration under Grant 
NAG5-9227 issued through the Office of Space Science, 
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and by National Science Foundation Grant AST9808980. lar secondary using adaptive optics observations with the 
We thank Sebastian Els for his attempt to image the stel- 4. 2m- William Herschel telescope. 
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Table 2 

McDonald Observatory Phase II radial velocities 



JD-2400000.0 


Radial Velocity 


a 




r — 1 1 

[ms L ] 


r — 1 1 

[ms L ] 


481 77 8696 

4:01 i i .OU^U 


3^ 39 

OO .OiJ 


33 09 


481 77 8773 

4:0 1 1 1 .O 1 IO 


90 QA 

LitJ . C/4: 


23.82 


48900 8^1 1 


97 79 


1 AO 


48994 7008 

4:OZZ4. 1 UUO 


^*3 Q*\ 

OO. c/O 


1 1 ^ 


489^0 fi090 


49 8^ 


1 3 no 


48484 R^Q^ 


*37 91 


1 zi nn 


48^94 8^04 

fOlJ^4:.OOW4: 


Q*3 1 


1 7 ^ 


48^^ 7Q^Q 

4:OUrJU. ( VOV 


10^1^ 


1 9 07 


48^^ 80^8 


1 1 01 


y.uo 


48fi07 7fi^8 


1 *38 1 8 
100. ±0 


1 8^ 


48fi44 FiAFi^ 


1 *37 70 


1 fi 41 


48894 0^.47 


904 ^0 


9D 9fi 


48894 QfiOQ 


1 Q9 79 
-iyz. ( z 


9zL 48 


488^9 Q7^ 


90Q 8*3 


9fl "?9 


488^ Q1 4Q 


91 9 fi4 


1 8 7fi 
10. |U 


48889 897*3 


90^ 97 


1 9 Qfi 


48Q01 7Q1 9 

■4:0 JU 1 . ( C/ _L Z 


99^ 14 

-ZiZO . J_4 


1 n 79 


48Q09 77fi4 


9*3^ 94 

-ZOu.Z4 


1 8 78 
10.10 


48Q4^ 744*3 

4:0C>40. ( 44:0 


94fi *38 


1 ^ Q1 


48Q71 fi^71 


9*38 fil 

ZiOO.Ul 


1 8 9"? 


48Q7^ fil 8Q 

4:0 C/ 1 O.U J-0*J 


90^ 78 


1 q ^.n 


4Q090 fi401 

4:C/UiiU.U4Ul 


91 8 Q7 


1 79 


4Q990 0,7*39 


1 (\7 9Q 
- id 1 .zy 


11 1Q 


409^8 8484 

4:>7i(J0.0404 


1 *39 8Q 


1 7 98 


4Q98fi 7QO*3 


181 fi^ 

101 -UO 


1 ^ Q9 


4Q*3^9 fi804 

4: c/ O O Z . U O U 4 


1 1 Q QQ 




40*380 ^871 

4:C/OOU.rJO ( _L 


Q7 9fi 


1 Q 91 


4:C/4UU.rJUrJU 


-1 1 ft 


17 10 


40^87 884^ 

4:c/UO ( .004rJ 


UO . (JO 


1 1 


40^87 8Q07 
4yjo * .oyu ( 


OA. .UO 


1 7 


4Qfi1 ^ 88QQ 


79 07 


14 

14:. OO 


4QB47 81 98 

4:^1)4 1 .Ol^O 


7^ ^^i 


1 00 


4QB70 71 1 8 


9S 9^i 

ZiO . Zi<J 


R 1 8 
0.10 


4Q70*3 fi78*3 

4:C/ 1 UO.U ( OO 


S1 

00. 1 


1 9 4 

J_ Zi .4: 


407*34 ^877 

4:c/ ( 04. rJO ( ( 


98 ^fi 
ZjO . 00 


1 ^ 9fi 
1 . zu 


4Q7BQ ^fi*3Q 


99 90 


1 1 4fi 

± ± .4U 




87 99 


1 9 9^ 

± Zi . LiO 


49946.9247 


103.39 


36.52 


49963.8567 


135.27 


20.53 


49993.8304 


127.05 


16.89 


50093.6243 


188.90 


8.42 


50124.6006 


222.54 


8.77 


50292.8908 


351.62 


13.65 


50354.8007 


389.66 


9.85 


50409.7363 


394.82 


14.13 


50480.6530 


420.67 


15.16 


50700.8995 


556.24 


16.33 


50768.8008 


575.22 


14.58 


50834.6049 


650.27 


12.11 



Hatzes et al. 



Table 3 

McDonald Observatory Phase III radial velocities 
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Table 4 

Velocity Offsets for the Datasets 



Dataset 


Velocity Offset 




[ms- 1 ] 


CFHT 


1294.6 ± 108 


McD Phase I 


2232.3 ± 108 


McD Phase II 


2040.8 ± 110 


McD Phase III 


864.7 ± 108 
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Table 5 

Binary Orbital Elements for 7 Cep 



Element This work Griffin ct al. (2002) 

Period (days) 20750.6579 ± 1568.6 24135 ± 349 

T (JD) 248429.03 ± 27.0 248625 ± 210 

Eccentricity 0.361 ± 0.023 0.389 ± 0.017 

uj (deg) 158.76 ± 1.2 166 ± 7 

Kl (kms -1 ) 1.82 ± 0.049 2.04 ± 0.10 

/(m) (solar masses) 0.0106 ± 0.0012 0.0166 ± 0.0025 

Semi-major axis (AU) 18.5 ± 1.1 20.3 ± 0.7 
Reduced x 2 (without planet) 4.36 



Table 6 

Orbital Elements for Planet around 7 Cep 



Element 

Period (days) 

T (JD) 
Eccentricity 

w (deg) 
Kl (ms- 1 ) 
/(to) (solar masses) 
Semi-major axis (AU) 
Reduced \ 2 
(Jcfht (ms -1 ) 

OPhasel (mS _1 ) 
OPhasell (mS -1 ) 
<?PhaseIII (mS -1 ) 



Value 



905.574 ± 3.08 

253121.925 ± 66.9 

0.12 ± 0.05 

49.6 ± 25.6 

27.50 ± 1.5 

(1.90 ± 0.3) xlO -9 

2.13 ± 0.05 

1.47 

15.3 

17.4 

15.8 

8.2 



Table 7 

Stellar Parameters (Fuhrmann 2003) 



Parameter 


Value 


T eff (days) 


4888 K 


log g 


3.33 


[Fe/H] 


+0.18 


M bol 


2.14 


Mass 


1.59 Mq 


Radius 


4.66 R G 


dnip 


13.79 pes 


dsp 


13.39 pes 
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Fig. 1. — The combined (planet + stellar) orbital solution to all RV data sets for 7 Cep. Circles represent the CFHT measurements 
of Walker et al. (1992). All other symbols are for data taken at McDonald Observatory. Triangles represent measurements using telluric 
O2 as the reference (Phase I), squares are for measurements using an iodine absorption cell (Phase II), and inverted triangles are also I2 
measurements, but with the large wavelength coverage 2dcoudc spectrograph (Phase III). 
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Frequency (c/d) 

Fig. 2. — The Lomb-Scargle periodogram of the combined RV measurements for 7 Cep after removal of the RV variations due to the 
stellar companion. (Top) The CFHT data alone. (Middle) The periodogram of the combined CHFT + McDonald Phases I, II data. (Bottom) 
Periodogram of all measurements. The false alarm probability (FAP) for the peak in each periodogram is shown in each panel. This was 
computed using the equation given in Scargle (1982). 
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Fig. 3. — The orbital solution for the planet (line) and the residual velocity measurements of the 4 data sets after subtracting the contribution 
due to the binary companion (points). 
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Fig. 4. — The phased residual RV measurements of CFHT, and McDonald Phase I III compared to the planet orbital solution (line). 
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Fig. 5. — The mean bisector span measurements (top) and the mean bisector curvature (bottom) measurements for 7 Cep phased to the 
planet orbital period. The dotted line marks zero value and the dashed line represents the extreme values of the radial velocity variations due 
to the planetary companion. 
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Fig. 6. — The Hipparcos photometry for 7 Cep over 1989.9 — 1993.17 phased to the planet orbital period. The crosses represent the 
individual measurements and the solid points phased-binned averages. The error bars indicate the rms scatter of values used in computing 
the binned average. 
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Fig. 7. — The Lomb-Scargle periodogram of the Hipparcos (daily averages) photometry. The vertical line marks the orbital frequency of 
the planet. 
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Fig. 9. — The S-index measurements from the McDonald Phase III data phased to the orbital period of the planet. Crosses represent 
the individual measurements, solid points phased-binned averages. The error bars represent the rms scatter of the data used for the binned 
averages. 
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Fig. 10. — The correlation between the RV and the Ca II A8662 A equivalent width (top) and the RV and the S-index measurements. A 
typical error bar is shown in the lower right of the figure. 
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Fig. 11. — The correlation between the RV and the McDonald Ca II S-index measurements measurements. 
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Fig. 13. — (Top) The Lomb-Scargle periodogram for the CHFT Ca II A8662 A measurements over the time span 1981 - 1986.5. (Bottom) 
The same for the time span 1986.5 - 1992. The vertical line marks the orbital frequency of the planet. 
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Fig. 14. — The 1986.5 — 1992 CFHT Ca II measurements phased to the period found in the periodogram analysis (top) and the planet 
period (middle). The bottom panel are the CFHT Ca II equivalent width variations from 1981-1986.5 phased to the planet period. 
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Fig. 15. — Lomb-Scargle periodogram of the RV-rcsiduals of 7 Cep after subtracting the binary and the planetary signal. The horizontal 
dashed lines show confidence levels with 10% and 1% falsc-alarm-probability (FAP). Obviously, no additional periodic signal above the noise 
level is present in the data. 



